3.885 \(\int \frac{x^2}{\left (a+b x^2+c x^4\right )^3} \, dx\)

Optimal. Leaf size=311 \[ -\frac{x \left (b+2 c x^2\right )}{4 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )^2}+\frac{x \left (c x^2 \left (20 a c+b^2\right )+b \left (8 a c+b^2\right )\right )}{8 a \left (b^2-4 a c\right )^2 \left (a+b x^2+c x^4\right )}+\frac{\sqrt{c} \left (\frac{b \left (b^2-52 a c\right )}{\sqrt{b^2-4 a c}}+20 a c+b^2\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{8 \sqrt{2} a \left (b^2-4 a c\right )^2 \sqrt{b-\sqrt{b^2-4 a c}}}+\frac{\sqrt{c} \left (-\frac{b \left (b^2-52 a c\right )}{\sqrt{b^2-4 a c}}+20 a c+b^2\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{\sqrt{b^2-4 a c}+b}}\right )}{8 \sqrt{2} a \left (b^2-4 a c\right )^2 \sqrt{\sqrt{b^2-4 a c}+b}} \]

[Out]

-(x*(b + 2*c*x^2))/(4*(b^2 - 4*a*c)*(a + b*x^2 + c*x^4)^2) + (x*(b*(b^2 + 8*a*c)
+ c*(b^2 + 20*a*c)*x^2))/(8*a*(b^2 - 4*a*c)^2*(a + b*x^2 + c*x^4)) + (Sqrt[c]*(
b^2 + 20*a*c + (b*(b^2 - 52*a*c))/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/
Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(8*Sqrt[2]*a*(b^2 - 4*a*c)^2*Sqrt[b - Sqrt[b^2 - 4
*a*c]]) + (Sqrt[c]*(b^2 + 20*a*c - (b*(b^2 - 52*a*c))/Sqrt[b^2 - 4*a*c])*ArcTan[
(Sqrt[2]*Sqrt[c]*x)/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(8*Sqrt[2]*a*(b^2 - 4*a*c)^2*S
qrt[b + Sqrt[b^2 - 4*a*c]])

_______________________________________________________________________________________

Rubi [A]  time = 1.50605, antiderivative size = 311, normalized size of antiderivative = 1., number of steps used = 5, number of rules used = 4, integrand size = 18, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.222 \[ -\frac{x \left (b+2 c x^2\right )}{4 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )^2}+\frac{x \left (c x^2 \left (20 a c+b^2\right )+b \left (8 a c+b^2\right )\right )}{8 a \left (b^2-4 a c\right )^2 \left (a+b x^2+c x^4\right )}+\frac{\sqrt{c} \left (\frac{b \left (b^2-52 a c\right )}{\sqrt{b^2-4 a c}}+20 a c+b^2\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{8 \sqrt{2} a \left (b^2-4 a c\right )^2 \sqrt{b-\sqrt{b^2-4 a c}}}+\frac{\sqrt{c} \left (-\frac{b \left (b^2-52 a c\right )}{\sqrt{b^2-4 a c}}+20 a c+b^2\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{\sqrt{b^2-4 a c}+b}}\right )}{8 \sqrt{2} a \left (b^2-4 a c\right )^2 \sqrt{\sqrt{b^2-4 a c}+b}} \]

Antiderivative was successfully verified.

[In]  Int[x^2/(a + b*x^2 + c*x^4)^3,x]

[Out]

-(x*(b + 2*c*x^2))/(4*(b^2 - 4*a*c)*(a + b*x^2 + c*x^4)^2) + (x*(b*(b^2 + 8*a*c)
+ c*(b^2 + 20*a*c)*x^2))/(8*a*(b^2 - 4*a*c)^2*(a + b*x^2 + c*x^4)) + (Sqrt[c]*(
b^2 + 20*a*c + (b*(b^2 - 52*a*c))/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/
Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(8*Sqrt[2]*a*(b^2 - 4*a*c)^2*Sqrt[b - Sqrt[b^2 - 4
*a*c]]) + (Sqrt[c]*(b^2 + 20*a*c - (b*(b^2 - 52*a*c))/Sqrt[b^2 - 4*a*c])*ArcTan[
(Sqrt[2]*Sqrt[c]*x)/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(8*Sqrt[2]*a*(b^2 - 4*a*c)^2*S
qrt[b + Sqrt[b^2 - 4*a*c]])

_______________________________________________________________________________________

Rubi in Sympy [A]  time = 81.3171, size = 291, normalized size = 0.94 \[ - \frac{x \left (b + 2 c x^{2}\right )}{4 \left (- 4 a c + b^{2}\right ) \left (a + b x^{2} + c x^{4}\right )^{2}} - \frac{\sqrt{2} \sqrt{c} \left (b \left (- 52 a c + b^{2}\right ) - \sqrt{- 4 a c + b^{2}} \left (20 a c + b^{2}\right )\right ) \operatorname{atan}{\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b + \sqrt{- 4 a c + b^{2}}}} \right )}}{16 a \sqrt{b + \sqrt{- 4 a c + b^{2}}} \left (- 4 a c + b^{2}\right )^{\frac{5}{2}}} + \frac{\sqrt{2} \sqrt{c} \left (b \left (- 52 a c + b^{2}\right ) + \sqrt{- 4 a c + b^{2}} \left (20 a c + b^{2}\right )\right ) \operatorname{atan}{\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b - \sqrt{- 4 a c + b^{2}}}} \right )}}{16 a \sqrt{b - \sqrt{- 4 a c + b^{2}}} \left (- 4 a c + b^{2}\right )^{\frac{5}{2}}} + \frac{x \left (b \left (8 a c + b^{2}\right ) + c x^{2} \left (20 a c + b^{2}\right )\right )}{8 a \left (- 4 a c + b^{2}\right )^{2} \left (a + b x^{2} + c x^{4}\right )} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  rubi_integrate(x**2/(c*x**4+b*x**2+a)**3,x)

[Out]

-x*(b + 2*c*x**2)/(4*(-4*a*c + b**2)*(a + b*x**2 + c*x**4)**2) - sqrt(2)*sqrt(c)
*(b*(-52*a*c + b**2) - sqrt(-4*a*c + b**2)*(20*a*c + b**2))*atan(sqrt(2)*sqrt(c)
*x/sqrt(b + sqrt(-4*a*c + b**2)))/(16*a*sqrt(b + sqrt(-4*a*c + b**2))*(-4*a*c +
b**2)**(5/2)) + sqrt(2)*sqrt(c)*(b*(-52*a*c + b**2) + sqrt(-4*a*c + b**2)*(20*a*
c + b**2))*atan(sqrt(2)*sqrt(c)*x/sqrt(b - sqrt(-4*a*c + b**2)))/(16*a*sqrt(b -
sqrt(-4*a*c + b**2))*(-4*a*c + b**2)**(5/2)) + x*(b*(8*a*c + b**2) + c*x**2*(20*
a*c + b**2))/(8*a*(-4*a*c + b**2)**2*(a + b*x**2 + c*x**4))

_______________________________________________________________________________________

Mathematica [A]  time = 1.54539, size = 334, normalized size = 1.07 \[ \frac{1}{16} \left (-\frac{4 x \left (b+2 c x^2\right )}{\left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )^2}+\frac{2 x \left (8 a b c+20 a c^2 x^2+b^3+b^2 c x^2\right )}{a \left (b^2-4 a c\right )^2 \left (a+b x^2+c x^4\right )}+\frac{\sqrt{2} \sqrt{c} \left (b^2 \sqrt{b^2-4 a c}+20 a c \sqrt{b^2-4 a c}-52 a b c+b^3\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{a \left (b^2-4 a c\right )^{5/2} \sqrt{b-\sqrt{b^2-4 a c}}}+\frac{\sqrt{2} \sqrt{c} \left (b^2 \sqrt{b^2-4 a c}+20 a c \sqrt{b^2-4 a c}+52 a b c-b^3\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{\sqrt{b^2-4 a c}+b}}\right )}{a \left (b^2-4 a c\right )^{5/2} \sqrt{\sqrt{b^2-4 a c}+b}}\right ) \]

Antiderivative was successfully verified.

[In]  Integrate[x^2/(a + b*x^2 + c*x^4)^3,x]

[Out]

((-4*x*(b + 2*c*x^2))/((b^2 - 4*a*c)*(a + b*x^2 + c*x^4)^2) + (2*x*(b^3 + 8*a*b*
c + b^2*c*x^2 + 20*a*c^2*x^2))/(a*(b^2 - 4*a*c)^2*(a + b*x^2 + c*x^4)) + (Sqrt[2
]*Sqrt[c]*(b^3 - 52*a*b*c + b^2*Sqrt[b^2 - 4*a*c] + 20*a*c*Sqrt[b^2 - 4*a*c])*Ar
cTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(a*(b^2 - 4*a*c)^(5/2)*Sq
rt[b - Sqrt[b^2 - 4*a*c]]) + (Sqrt[2]*Sqrt[c]*(-b^3 + 52*a*b*c + b^2*Sqrt[b^2 -
4*a*c] + 20*a*c*Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b + Sqrt[b^2
- 4*a*c]]])/(a*(b^2 - 4*a*c)^(5/2)*Sqrt[b + Sqrt[b^2 - 4*a*c]]))/16

_______________________________________________________________________________________

Maple [B]  time = 0.165, size = 2958, normalized size = 9.5 \[ \text{output too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  int(x^2/(c*x^4+b*x^2+a)^3,x)

[Out]

-9*c^3/(-4*a*c+b^2)^2/(4*a*c-b^2)^2*a*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*a
rctan(c*x*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b^2-20*c^4/(-4*a*c+b^2)^2/(4
*a*c-b^2)^2*a^2*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctanh(c*x*2^(1/2)/((
-b+(-4*a*c+b^2)^(1/2))*c)^(1/2))-4*c^3/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2/(x^2+1/2
*b/c-1/2/c*(-4*a*c+b^2)^(1/2))^2*a^2*x^3*b+3*c^2/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^
2/(x^2+1/2*b/c-1/2/c*(-4*a*c+b^2)^(1/2))^2*a*b^3*x^3-42*c^2/(-4*a*c+b^2)^(5/2)/(
4*a*c-b^2)^2/(x^2+1/2*b/c-1/2/c*(-4*a*c+b^2)^(1/2))^2*a^2*b^2*x+21/2*c/(-4*a*c+b
^2)^(5/2)/(4*a*c-b^2)^2/(x^2+1/2*b/c-1/2/c*(-4*a*c+b^2)^(1/2))^2*x*a*b^4-3/4*c^2
/(-4*a*c+b^2)^2/(4*a*c-b^2)^2*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctanh(
c*x*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b^4+20*c^4/(-4*a*c+b^2)^2/(4*a*c-
b^2)^2*a^2*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2)/((b+(-4*a
*c+b^2)^(1/2))*c)^(1/2))-9*c^2/(-4*a*c+b^2)^2/(4*a*c-b^2)^2/(x^2+1/2/c*(-4*a*c+b
^2)^(1/2)+1/2*b/c)^2*a*x^3*b^2+9*c^3/(-4*a*c+b^2)^2/(4*a*c-b^2)^2*a*2^(1/2)/((-b
+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctanh(c*x*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1
/2))*b^2+3/4/(-4*a*c+b^2)^2/(4*a*c-b^2)^2/(x^2+1/2*b/c-1/2/c*(-4*a*c+b^2)^(1/2))
^2*x*b^5-7/8/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2/(x^2+1/2*b/c-1/2/c*(-4*a*c+b^2)^(1
/2))^2*x*b^6+7/8/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2/(x^2+1/2/c*(-4*a*c+b^2)^(1/2)+
1/2*b/c)^2*x*b^6+12*c^2/(-4*a*c+b^2)^2/(4*a*c-b^2)^2/(x^2+1/2*b/c-1/2/c*(-4*a*c+
b^2)^(1/2))^2*x*a^2*b-6*c/(-4*a*c+b^2)^2/(4*a*c-b^2)^2/(x^2+1/2*b/c-1/2/c*(-4*a*
c+b^2)^(1/2))^2*x*a*b^3-9*c^2/(-4*a*c+b^2)^2/(4*a*c-b^2)^2/(x^2+1/2*b/c-1/2/c*(-
4*a*c+b^2)^(1/2))^2*a*x^3*b^2-3*c^2/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2/(x^2+1/2/c*
(-4*a*c+b^2)^(1/2)+1/2*b/c)^2*a*b^3*x^3+42*c^2/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2/
(x^2+1/2/c*(-4*a*c+b^2)^(1/2)+1/2*b/c)^2*a^2*b^2*x-21/2*c/(-4*a*c+b^2)^(5/2)/(4*
a*c-b^2)^2/(x^2+1/2/c*(-4*a*c+b^2)^(1/2)+1/2*b/c)^2*x*a*b^4+4*c^3/(-4*a*c+b^2)^(
5/2)/(4*a*c-b^2)^2/(x^2+1/2/c*(-4*a*c+b^2)^(1/2)+1/2*b/c)^2*a^2*x^3*b+3/4/(-4*a*
c+b^2)^2/(4*a*c-b^2)^2/(x^2+1/2/c*(-4*a*c+b^2)^(1/2)+1/2*b/c)^2*x*b^5+15/4*c^2/(
-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(
c*x*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b^5+15/4*c^2/(-4*a*c+b^2)^(5/2)/(4
*a*c-b^2)^2*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctanh(c*x*2^(1/2)/((-b+(
-4*a*c+b^2)^(1/2))*c)^(1/2))*b^5+3/4*c^2/(-4*a*c+b^2)^2/(4*a*c-b^2)^2*2^(1/2)/((
b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/
2))*b^4+12*c^2/(-4*a*c+b^2)^2/(4*a*c-b^2)^2/(x^2+1/2/c*(-4*a*c+b^2)^(1/2)+1/2*b/
c)^2*x*a^2*b-6*c/(-4*a*c+b^2)^2/(4*a*c-b^2)^2/(x^2+1/2/c*(-4*a*c+b^2)^(1/2)+1/2*
b/c)^2*x*a*b^3+52*c^4/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2*a^2*2^(1/2)/((b+(-4*a*c+b
^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b-27*c^
3/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2*a*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*ar
ctan(c*x*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b^3-1/16*c/(-4*a*c+b^2)^(5/2)
/(4*a*c-b^2)^2/a*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2)/((b
+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b^7+52*c^4/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2*a^2*2
^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctanh(c*x*2^(1/2)/((-b+(-4*a*c+b^2)^(
1/2))*c)^(1/2))*b-27*c^3/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2*a*2^(1/2)/((-b+(-4*a*c
+b^2)^(1/2))*c)^(1/2)*arctanh(c*x*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b^3
-1/16*c/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2/a*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(
1/2)*arctanh(c*x*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b^7+1/16*c/(-4*a*c+b
^2)^2/(4*a*c-b^2)^2/a*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2
)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b^6-1/16*c/(-4*a*c+b^2)^2/(4*a*c-b^2)^2/a*2^
(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctanh(c*x*2^(1/2)/((-b+(-4*a*c+b^2)^(1
/2))*c)^(1/2))*b^6+3/4*c/(-4*a*c+b^2)^2/(4*a*c-b^2)^2/(x^2+1/2/c*(-4*a*c+b^2)^(1
/2)+1/2*b/c)^2*x^3*b^4+1/16/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2/(x^2+1/2*b/c-1/2/c*
(-4*a*c+b^2)^(1/2))^2/a*x^3*b^7+20*c^3/(-4*a*c+b^2)^2/(4*a*c-b^2)^2/(x^2+1/2/c*(
-4*a*c+b^2)^(1/2)+1/2*b/c)^2*a^2*x^3+3/4*c/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2/(x^2
+1/2/c*(-4*a*c+b^2)^(1/2)+1/2*b/c)^2*x^3*b^5-56*c^3/(-4*a*c+b^2)^(5/2)/(4*a*c-b^
2)^2/(x^2+1/2/c*(-4*a*c+b^2)^(1/2)+1/2*b/c)^2*x*a^3+3/4*c/(-4*a*c+b^2)^2/(4*a*c-
b^2)^2/(x^2+1/2*b/c-1/2/c*(-4*a*c+b^2)^(1/2))^2*x^3*b^4+20*c^3/(-4*a*c+b^2)^2/(4
*a*c-b^2)^2/(x^2+1/2*b/c-1/2/c*(-4*a*c+b^2)^(1/2))^2*a^2*x^3-3/4*c/(-4*a*c+b^2)^
(5/2)/(4*a*c-b^2)^2/(x^2+1/2*b/c-1/2/c*(-4*a*c+b^2)^(1/2))^2*x^3*b^5+56*c^3/(-4*
a*c+b^2)^(5/2)/(4*a*c-b^2)^2/(x^2+1/2*b/c-1/2/c*(-4*a*c+b^2)^(1/2))^2*x*a^3+1/16
/(-4*a*c+b^2)^2/(4*a*c-b^2)^2/(x^2+1/2/c*(-4*a*c+b^2)^(1/2)+1/2*b/c)^2/a*x^3*b^6
+1/16/(-4*a*c+b^2)^2/(4*a*c-b^2)^2/(x^2+1/2*b/c-1/2/c*(-4*a*c+b^2)^(1/2))^2/a*x^
3*b^6-1/16/(-4*a*c+b^2)^(5/2)/(4*a*c-b^2)^2/(x^2+1/2/c*(-4*a*c+b^2)^(1/2)+1/2*b/
c)^2/a*x^3*b^7

_______________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \[ \frac{{\left (b^{2} c^{2} + 20 \, a c^{3}\right )} x^{7} + 2 \,{\left (b^{3} c + 14 \, a b c^{2}\right )} x^{5} +{\left (b^{4} + 5 \, a b^{2} c + 36 \, a^{2} c^{2}\right )} x^{3} -{\left (a b^{3} - 16 \, a^{2} b c\right )} x}{8 \,{\left ({\left (a b^{4} c^{2} - 8 \, a^{2} b^{2} c^{3} + 16 \, a^{3} c^{4}\right )} x^{8} + a^{3} b^{4} - 8 \, a^{4} b^{2} c + 16 \, a^{5} c^{2} + 2 \,{\left (a b^{5} c - 8 \, a^{2} b^{3} c^{2} + 16 \, a^{3} b c^{3}\right )} x^{6} +{\left (a b^{6} - 6 \, a^{2} b^{4} c + 32 \, a^{4} c^{3}\right )} x^{4} + 2 \,{\left (a^{2} b^{5} - 8 \, a^{3} b^{3} c + 16 \, a^{4} b c^{2}\right )} x^{2}\right )}} + \frac{\int \frac{b^{3} - 16 \, a b c +{\left (b^{2} c + 20 \, a c^{2}\right )} x^{2}}{c x^{4} + b x^{2} + a}\,{d x}}{8 \,{\left (a b^{4} - 8 \, a^{2} b^{2} c + 16 \, a^{3} c^{2}\right )}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(x^2/(c*x^4 + b*x^2 + a)^3,x, algorithm="maxima")

[Out]

1/8*((b^2*c^2 + 20*a*c^3)*x^7 + 2*(b^3*c + 14*a*b*c^2)*x^5 + (b^4 + 5*a*b^2*c +
36*a^2*c^2)*x^3 - (a*b^3 - 16*a^2*b*c)*x)/((a*b^4*c^2 - 8*a^2*b^2*c^3 + 16*a^3*c
^4)*x^8 + a^3*b^4 - 8*a^4*b^2*c + 16*a^5*c^2 + 2*(a*b^5*c - 8*a^2*b^3*c^2 + 16*a
^3*b*c^3)*x^6 + (a*b^6 - 6*a^2*b^4*c + 32*a^4*c^3)*x^4 + 2*(a^2*b^5 - 8*a^3*b^3*
c + 16*a^4*b*c^2)*x^2) + 1/8*integrate((b^3 - 16*a*b*c + (b^2*c + 20*a*c^2)*x^2)
/(c*x^4 + b*x^2 + a), x)/(a*b^4 - 8*a^2*b^2*c + 16*a^3*c^2)

_______________________________________________________________________________________

Fricas [A]  time = 0.42304, size = 5099, normalized size = 16.4 \[ \text{result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(x^2/(c*x^4 + b*x^2 + a)^3,x, algorithm="fricas")

[Out]

1/16*(2*(b^2*c^2 + 20*a*c^3)*x^7 + 4*(b^3*c + 14*a*b*c^2)*x^5 + 2*(b^4 + 5*a*b^2
*c + 36*a^2*c^2)*x^3 + sqrt(1/2)*((a*b^4*c^2 - 8*a^2*b^2*c^3 + 16*a^3*c^4)*x^8 +
a^3*b^4 - 8*a^4*b^2*c + 16*a^5*c^2 + 2*(a*b^5*c - 8*a^2*b^3*c^2 + 16*a^3*b*c^3)
*x^6 + (a*b^6 - 6*a^2*b^4*c + 32*a^4*c^3)*x^4 + 2*(a^2*b^5 - 8*a^3*b^3*c + 16*a^
4*b*c^2)*x^2)*sqrt(-(b^7 - 35*a*b^5*c + 280*a^2*b^3*c^2 + 1680*a^3*b*c^3 + (a^3*
b^10 - 20*a^4*b^8*c + 160*a^5*b^6*c^2 - 640*a^6*b^4*c^3 + 1280*a^7*b^2*c^4 - 102
4*a^8*c^5)*sqrt((b^4 - 50*a*b^2*c + 625*a^2*c^2)/(a^6*b^10 - 20*a^7*b^8*c + 160*
a^8*b^6*c^2 - 640*a^9*b^4*c^3 + 1280*a^10*b^2*c^4 - 1024*a^11*c^5)))/(a^3*b^10 -
20*a^4*b^8*c + 160*a^5*b^6*c^2 - 640*a^6*b^4*c^3 + 1280*a^7*b^2*c^4 - 1024*a^8*
c^5))*log((35*b^6*c^2 - 1491*a*b^4*c^3 + 15000*a^2*b^2*c^4 + 10000*a^3*c^5)*x +
1/2*sqrt(1/2)*(b^11 - 53*a*b^9*c + 940*a^2*b^7*c^2 - 6832*a^3*b^5*c^3 + 21824*a^
4*b^3*c^4 - 25600*a^5*b*c^5 - (a^3*b^14 - 38*a^4*b^12*c + 480*a^5*b^10*c^2 - 272
0*a^6*b^8*c^3 + 6400*a^7*b^6*c^4 + 1536*a^8*b^4*c^5 - 32768*a^9*b^2*c^6 + 40960*
a^10*c^7)*sqrt((b^4 - 50*a*b^2*c + 625*a^2*c^2)/(a^6*b^10 - 20*a^7*b^8*c + 160*a
^8*b^6*c^2 - 640*a^9*b^4*c^3 + 1280*a^10*b^2*c^4 - 1024*a^11*c^5)))*sqrt(-(b^7 -
35*a*b^5*c + 280*a^2*b^3*c^2 + 1680*a^3*b*c^3 + (a^3*b^10 - 20*a^4*b^8*c + 160*
a^5*b^6*c^2 - 640*a^6*b^4*c^3 + 1280*a^7*b^2*c^4 - 1024*a^8*c^5)*sqrt((b^4 - 50*
a*b^2*c + 625*a^2*c^2)/(a^6*b^10 - 20*a^7*b^8*c + 160*a^8*b^6*c^2 - 640*a^9*b^4*
c^3 + 1280*a^10*b^2*c^4 - 1024*a^11*c^5)))/(a^3*b^10 - 20*a^4*b^8*c + 160*a^5*b^
6*c^2 - 640*a^6*b^4*c^3 + 1280*a^7*b^2*c^4 - 1024*a^8*c^5))) - sqrt(1/2)*((a*b^4
*c^2 - 8*a^2*b^2*c^3 + 16*a^3*c^4)*x^8 + a^3*b^4 - 8*a^4*b^2*c + 16*a^5*c^2 + 2*
(a*b^5*c - 8*a^2*b^3*c^2 + 16*a^3*b*c^3)*x^6 + (a*b^6 - 6*a^2*b^4*c + 32*a^4*c^3
)*x^4 + 2*(a^2*b^5 - 8*a^3*b^3*c + 16*a^4*b*c^2)*x^2)*sqrt(-(b^7 - 35*a*b^5*c +
280*a^2*b^3*c^2 + 1680*a^3*b*c^3 + (a^3*b^10 - 20*a^4*b^8*c + 160*a^5*b^6*c^2 -
640*a^6*b^4*c^3 + 1280*a^7*b^2*c^4 - 1024*a^8*c^5)*sqrt((b^4 - 50*a*b^2*c + 625*
a^2*c^2)/(a^6*b^10 - 20*a^7*b^8*c + 160*a^8*b^6*c^2 - 640*a^9*b^4*c^3 + 1280*a^1
0*b^2*c^4 - 1024*a^11*c^5)))/(a^3*b^10 - 20*a^4*b^8*c + 160*a^5*b^6*c^2 - 640*a^
6*b^4*c^3 + 1280*a^7*b^2*c^4 - 1024*a^8*c^5))*log((35*b^6*c^2 - 1491*a*b^4*c^3 +
15000*a^2*b^2*c^4 + 10000*a^3*c^5)*x - 1/2*sqrt(1/2)*(b^11 - 53*a*b^9*c + 940*a
^2*b^7*c^2 - 6832*a^3*b^5*c^3 + 21824*a^4*b^3*c^4 - 25600*a^5*b*c^5 - (a^3*b^14
- 38*a^4*b^12*c + 480*a^5*b^10*c^2 - 2720*a^6*b^8*c^3 + 6400*a^7*b^6*c^4 + 1536*
a^8*b^4*c^5 - 32768*a^9*b^2*c^6 + 40960*a^10*c^7)*sqrt((b^4 - 50*a*b^2*c + 625*a
^2*c^2)/(a^6*b^10 - 20*a^7*b^8*c + 160*a^8*b^6*c^2 - 640*a^9*b^4*c^3 + 1280*a^10
*b^2*c^4 - 1024*a^11*c^5)))*sqrt(-(b^7 - 35*a*b^5*c + 280*a^2*b^3*c^2 + 1680*a^3
*b*c^3 + (a^3*b^10 - 20*a^4*b^8*c + 160*a^5*b^6*c^2 - 640*a^6*b^4*c^3 + 1280*a^7
*b^2*c^4 - 1024*a^8*c^5)*sqrt((b^4 - 50*a*b^2*c + 625*a^2*c^2)/(a^6*b^10 - 20*a^
7*b^8*c + 160*a^8*b^6*c^2 - 640*a^9*b^4*c^3 + 1280*a^10*b^2*c^4 - 1024*a^11*c^5)
))/(a^3*b^10 - 20*a^4*b^8*c + 160*a^5*b^6*c^2 - 640*a^6*b^4*c^3 + 1280*a^7*b^2*c
^4 - 1024*a^8*c^5))) + sqrt(1/2)*((a*b^4*c^2 - 8*a^2*b^2*c^3 + 16*a^3*c^4)*x^8 +
a^3*b^4 - 8*a^4*b^2*c + 16*a^5*c^2 + 2*(a*b^5*c - 8*a^2*b^3*c^2 + 16*a^3*b*c^3)
*x^6 + (a*b^6 - 6*a^2*b^4*c + 32*a^4*c^3)*x^4 + 2*(a^2*b^5 - 8*a^3*b^3*c + 16*a^
4*b*c^2)*x^2)*sqrt(-(b^7 - 35*a*b^5*c + 280*a^2*b^3*c^2 + 1680*a^3*b*c^3 - (a^3*
b^10 - 20*a^4*b^8*c + 160*a^5*b^6*c^2 - 640*a^6*b^4*c^3 + 1280*a^7*b^2*c^4 - 102
4*a^8*c^5)*sqrt((b^4 - 50*a*b^2*c + 625*a^2*c^2)/(a^6*b^10 - 20*a^7*b^8*c + 160*
a^8*b^6*c^2 - 640*a^9*b^4*c^3 + 1280*a^10*b^2*c^4 - 1024*a^11*c^5)))/(a^3*b^10 -
20*a^4*b^8*c + 160*a^5*b^6*c^2 - 640*a^6*b^4*c^3 + 1280*a^7*b^2*c^4 - 1024*a^8*
c^5))*log((35*b^6*c^2 - 1491*a*b^4*c^3 + 15000*a^2*b^2*c^4 + 10000*a^3*c^5)*x +
1/2*sqrt(1/2)*(b^11 - 53*a*b^9*c + 940*a^2*b^7*c^2 - 6832*a^3*b^5*c^3 + 21824*a^
4*b^3*c^4 - 25600*a^5*b*c^5 + (a^3*b^14 - 38*a^4*b^12*c + 480*a^5*b^10*c^2 - 272
0*a^6*b^8*c^3 + 6400*a^7*b^6*c^4 + 1536*a^8*b^4*c^5 - 32768*a^9*b^2*c^6 + 40960*
a^10*c^7)*sqrt((b^4 - 50*a*b^2*c + 625*a^2*c^2)/(a^6*b^10 - 20*a^7*b^8*c + 160*a
^8*b^6*c^2 - 640*a^9*b^4*c^3 + 1280*a^10*b^2*c^4 - 1024*a^11*c^5)))*sqrt(-(b^7 -
35*a*b^5*c + 280*a^2*b^3*c^2 + 1680*a^3*b*c^3 - (a^3*b^10 - 20*a^4*b^8*c + 160*
a^5*b^6*c^2 - 640*a^6*b^4*c^3 + 1280*a^7*b^2*c^4 - 1024*a^8*c^5)*sqrt((b^4 - 50*
a*b^2*c + 625*a^2*c^2)/(a^6*b^10 - 20*a^7*b^8*c + 160*a^8*b^6*c^2 - 640*a^9*b^4*
c^3 + 1280*a^10*b^2*c^4 - 1024*a^11*c^5)))/(a^3*b^10 - 20*a^4*b^8*c + 160*a^5*b^
6*c^2 - 640*a^6*b^4*c^3 + 1280*a^7*b^2*c^4 - 1024*a^8*c^5))) - sqrt(1/2)*((a*b^4
*c^2 - 8*a^2*b^2*c^3 + 16*a^3*c^4)*x^8 + a^3*b^4 - 8*a^4*b^2*c + 16*a^5*c^2 + 2*
(a*b^5*c - 8*a^2*b^3*c^2 + 16*a^3*b*c^3)*x^6 + (a*b^6 - 6*a^2*b^4*c + 32*a^4*c^3
)*x^4 + 2*(a^2*b^5 - 8*a^3*b^3*c + 16*a^4*b*c^2)*x^2)*sqrt(-(b^7 - 35*a*b^5*c +
280*a^2*b^3*c^2 + 1680*a^3*b*c^3 - (a^3*b^10 - 20*a^4*b^8*c + 160*a^5*b^6*c^2 -
640*a^6*b^4*c^3 + 1280*a^7*b^2*c^4 - 1024*a^8*c^5)*sqrt((b^4 - 50*a*b^2*c + 625*
a^2*c^2)/(a^6*b^10 - 20*a^7*b^8*c + 160*a^8*b^6*c^2 - 640*a^9*b^4*c^3 + 1280*a^1
0*b^2*c^4 - 1024*a^11*c^5)))/(a^3*b^10 - 20*a^4*b^8*c + 160*a^5*b^6*c^2 - 640*a^
6*b^4*c^3 + 1280*a^7*b^2*c^4 - 1024*a^8*c^5))*log((35*b^6*c^2 - 1491*a*b^4*c^3 +
15000*a^2*b^2*c^4 + 10000*a^3*c^5)*x - 1/2*sqrt(1/2)*(b^11 - 53*a*b^9*c + 940*a
^2*b^7*c^2 - 6832*a^3*b^5*c^3 + 21824*a^4*b^3*c^4 - 25600*a^5*b*c^5 + (a^3*b^14
- 38*a^4*b^12*c + 480*a^5*b^10*c^2 - 2720*a^6*b^8*c^3 + 6400*a^7*b^6*c^4 + 1536*
a^8*b^4*c^5 - 32768*a^9*b^2*c^6 + 40960*a^10*c^7)*sqrt((b^4 - 50*a*b^2*c + 625*a
^2*c^2)/(a^6*b^10 - 20*a^7*b^8*c + 160*a^8*b^6*c^2 - 640*a^9*b^4*c^3 + 1280*a^10
*b^2*c^4 - 1024*a^11*c^5)))*sqrt(-(b^7 - 35*a*b^5*c + 280*a^2*b^3*c^2 + 1680*a^3
*b*c^3 - (a^3*b^10 - 20*a^4*b^8*c + 160*a^5*b^6*c^2 - 640*a^6*b^4*c^3 + 1280*a^7
*b^2*c^4 - 1024*a^8*c^5)*sqrt((b^4 - 50*a*b^2*c + 625*a^2*c^2)/(a^6*b^10 - 20*a^
7*b^8*c + 160*a^8*b^6*c^2 - 640*a^9*b^4*c^3 + 1280*a^10*b^2*c^4 - 1024*a^11*c^5)
))/(a^3*b^10 - 20*a^4*b^8*c + 160*a^5*b^6*c^2 - 640*a^6*b^4*c^3 + 1280*a^7*b^2*c
^4 - 1024*a^8*c^5))) - 2*(a*b^3 - 16*a^2*b*c)*x)/((a*b^4*c^2 - 8*a^2*b^2*c^3 + 1
6*a^3*c^4)*x^8 + a^3*b^4 - 8*a^4*b^2*c + 16*a^5*c^2 + 2*(a*b^5*c - 8*a^2*b^3*c^2
+ 16*a^3*b*c^3)*x^6 + (a*b^6 - 6*a^2*b^4*c + 32*a^4*c^3)*x^4 + 2*(a^2*b^5 - 8*a
^3*b^3*c + 16*a^4*b*c^2)*x^2)

_______________________________________________________________________________________

Sympy [A]  time = 53.3318, size = 733, normalized size = 2.36 \[ \frac{x^{7} \left (20 a c^{3} + b^{2} c^{2}\right ) + x^{5} \left (28 a b c^{2} + 2 b^{3} c\right ) + x^{3} \left (36 a^{2} c^{2} + 5 a b^{2} c + b^{4}\right ) + x \left (16 a^{2} b c - a b^{3}\right )}{128 a^{5} c^{2} - 64 a^{4} b^{2} c + 8 a^{3} b^{4} + x^{8} \left (128 a^{3} c^{4} - 64 a^{2} b^{2} c^{3} + 8 a b^{4} c^{2}\right ) + x^{6} \left (256 a^{3} b c^{3} - 128 a^{2} b^{3} c^{2} + 16 a b^{5} c\right ) + x^{4} \left (256 a^{4} c^{3} - 48 a^{2} b^{4} c + 8 a b^{6}\right ) + x^{2} \left (256 a^{4} b c^{2} - 128 a^{3} b^{3} c + 16 a^{2} b^{5}\right )} + \operatorname{RootSum}{\left (t^{4} \left (68719476736 a^{13} c^{10} - 171798691840 a^{12} b^{2} c^{9} + 193273528320 a^{11} b^{4} c^{8} - 128849018880 a^{10} b^{6} c^{7} + 56371445760 a^{9} b^{8} c^{6} - 16911433728 a^{8} b^{10} c^{5} + 3523215360 a^{7} b^{12} c^{4} - 503316480 a^{6} b^{14} c^{3} + 47185920 a^{5} b^{16} c^{2} - 2621440 a^{4} b^{18} c + 65536 a^{3} b^{20}\right ) + t^{2} \left (- 440401920 a^{8} b c^{8} + 477102080 a^{7} b^{3} c^{7} - 174325760 a^{6} b^{5} c^{6} + 11206656 a^{5} b^{7} c^{5} + 8929280 a^{4} b^{9} c^{4} - 2600960 a^{3} b^{11} c^{3} + 291840 a^{2} b^{13} c^{2} - 14080 a b^{15} c + 256 b^{17}\right ) + 160000 a^{4} c^{7} + 492800 a^{3} b^{2} c^{6} + 351456 a^{2} b^{4} c^{5} - 43120 a b^{6} c^{4} + 1225 b^{8} c^{3}, \left ( t \mapsto t \log{\left (x + \frac{167772160 t^{3} a^{10} c^{7} - 134217728 t^{3} a^{9} b^{2} c^{6} + 6291456 t^{3} a^{8} b^{4} c^{5} + 26214400 t^{3} a^{7} b^{6} c^{4} - 11141120 t^{3} a^{6} b^{8} c^{3} + 1966080 t^{3} a^{5} b^{10} c^{2} - 155648 t^{3} a^{4} b^{12} c + 4096 t^{3} a^{3} b^{14} - 742400 t a^{5} b c^{5} - 156928 t a^{4} b^{3} c^{4} - 70336 t a^{3} b^{5} c^{3} + 14480 t a^{2} b^{7} c^{2} - 848 t a b^{9} c + 16 t b^{11}}{10000 a^{3} c^{5} + 15000 a^{2} b^{2} c^{4} - 1491 a b^{4} c^{3} + 35 b^{6} c^{2}} \right )} \right )\right )} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(x**2/(c*x**4+b*x**2+a)**3,x)

[Out]

(x**7*(20*a*c**3 + b**2*c**2) + x**5*(28*a*b*c**2 + 2*b**3*c) + x**3*(36*a**2*c*
*2 + 5*a*b**2*c + b**4) + x*(16*a**2*b*c - a*b**3))/(128*a**5*c**2 - 64*a**4*b**
2*c + 8*a**3*b**4 + x**8*(128*a**3*c**4 - 64*a**2*b**2*c**3 + 8*a*b**4*c**2) + x
**6*(256*a**3*b*c**3 - 128*a**2*b**3*c**2 + 16*a*b**5*c) + x**4*(256*a**4*c**3 -
48*a**2*b**4*c + 8*a*b**6) + x**2*(256*a**4*b*c**2 - 128*a**3*b**3*c + 16*a**2*
b**5)) + RootSum(_t**4*(68719476736*a**13*c**10 - 171798691840*a**12*b**2*c**9 +
193273528320*a**11*b**4*c**8 - 128849018880*a**10*b**6*c**7 + 56371445760*a**9*
b**8*c**6 - 16911433728*a**8*b**10*c**5 + 3523215360*a**7*b**12*c**4 - 503316480
*a**6*b**14*c**3 + 47185920*a**5*b**16*c**2 - 2621440*a**4*b**18*c + 65536*a**3*
b**20) + _t**2*(-440401920*a**8*b*c**8 + 477102080*a**7*b**3*c**7 - 174325760*a*
*6*b**5*c**6 + 11206656*a**5*b**7*c**5 + 8929280*a**4*b**9*c**4 - 2600960*a**3*b
**11*c**3 + 291840*a**2*b**13*c**2 - 14080*a*b**15*c + 256*b**17) + 160000*a**4*
c**7 + 492800*a**3*b**2*c**6 + 351456*a**2*b**4*c**5 - 43120*a*b**6*c**4 + 1225*
b**8*c**3, Lambda(_t, _t*log(x + (167772160*_t**3*a**10*c**7 - 134217728*_t**3*a
**9*b**2*c**6 + 6291456*_t**3*a**8*b**4*c**5 + 26214400*_t**3*a**7*b**6*c**4 - 1
1141120*_t**3*a**6*b**8*c**3 + 1966080*_t**3*a**5*b**10*c**2 - 155648*_t**3*a**4
*b**12*c + 4096*_t**3*a**3*b**14 - 742400*_t*a**5*b*c**5 - 156928*_t*a**4*b**3*c
**4 - 70336*_t*a**3*b**5*c**3 + 14480*_t*a**2*b**7*c**2 - 848*_t*a*b**9*c + 16*_
t*b**11)/(10000*a**3*c**5 + 15000*a**2*b**2*c**4 - 1491*a*b**4*c**3 + 35*b**6*c*
*2))))

_______________________________________________________________________________________

GIAC/XCAS [F(-2)]  time = 0., size = 0, normalized size = 0. \[ \text{Exception raised: TypeError} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]  integrate(x^2/(c*x^4 + b*x^2 + a)^3,x, algorithm="giac")

[Out]

Exception raised: TypeError